Interaction of Porosity with a Planar Solid/Liquid Interface 
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In this article, an investigation of the interaction between gas porosity and a planar solid/liquid (SL) 
interface is reported. A two-dimensional numerical model able to accurately track sharp SL interfaces 
during solidification of pure metals and alloys is proposed. The finite-difference method and a rectangular 
undeformed grid are used for computation. The SL interface is described through the points of inter- 
section with the grid lines. Its motion is determined by the thermal and solute gradients at each par- 
ticular point. Changes of the interface temperature because of capillarity or solute redistribution as well 
as any perturbation of the thermal and solute field produced by the presence of non-metallic inclu- 
sions can be computed. To validate the model, the dynamics of the interaction between a gas pore and 
a solidification front in metal alloys was observed using a state of the art X-ray transmission micro- 
scope (XTM). The experiments included observation of the distortion of the SL interface near a pore, 
real-time measurements of the growth rate, and the change in shape of the porosity during interaction 
with the SL interface in pure A1 and Al-0.25 wt pet Au alloy. In addition, porosity-induced solute seg- 
regation patterns surrounding a pore were also quantified. 


I. INTRODUCTION 

GAS evolution during solidification is responsible for 
two casting defects: macroporosity (gas porosity ) and micro- 
shrinkage (microporosity). While both defects have a signifi- 
cant influence on mechanical properties, their mechanism of 
formation is quite different. 1 1] Macroporosity results when gas 
is rejected from the liquid and is entrapped in the solidifying 
metals as spherical gas bubbles of millimeter size. Micro- 
shrinkage occurs when liquid metal cannot reach interden- 
dritic areas during casting solidification, and is caused by a 
combination of shrinkage and gas evolution. It is a stand- 
ard foundry defect for mushy-freezing alloys, such as alu- 
minum alloys, steel, superalloys, brass, bronze, and cast iron. 
Its size is of the order of 10 to 100 pm. 

Previous experimental studies of porosity evolution in 
metals have primarily relied on quenching of directionally 
solidified samples followed by postsolidification analysis of 
the combined effect of gas and shrinkage porosity. For exam- 
ple, experiments performed on Cr, Ni, and Mo alloyed steel 121 
revealed that elongated or ellipsoidal macroporosity occurred 
only when a nucleation agent such as A1 2 0 3 was added to 
the melt. In the absence of nucleants, irrespective of the 
hydrogen content, only microshrinkage was obtained. The 
effect of gravity on pore distribution was studied by Kim 
et alP ] in a A356 Al-Si alloy. They demonstrated that for 
upward solidification, the distribution of pores was con- 
centrated at the hot (upper) end of the sample, implying that 
terminal floatation velocity was the dominant factor. Floata- 
tion did not permit interaction between the pores and the 
solid/liquid (SL) interface. During downward solidification, 
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both macro- and microshrinkage were found entrapped in 
the mushy zone. However, postsolidification analysis is inca- 
pable of either capturing the dynamics of porosity evolution 
or their interaction with the solidification front. 

To overcome this lack of real-time data, Lee and Hunt- 41 
used a microfocus X-ray system with a resolution of 25 pm 
to study the real-time interaction between a pore and primary 
dendrites in Al-10 wt pet Cu samples. Their experimental 
observation revealed that the shape of the pore undergoes a 
transition from spherical to ellipsoidal once the pore is impinged 
upon by the dendritic interface. With an increase in growth 
velocity, the primary dendrite arm spacing decreases and fur- 
ther elongates the pore in the growth direction. However, this 
investigation did not address the issue of solutal segregation 
that may be induced by the porosity during solidification. 

Real-time observations on the interaction between the solidi- 
fication front and a pore have also been performed on trans- 
parent organic (analogue) materials. Jamgotchian et al. l5] 
demonstrated that in the succinonitrile-acetone system, gas 
bubbles that formed at the planar SL interface grew in a cou- 
pled fashion into elongated shapes called duplexes. Experi- 
ments with the faceted salol system have shown that there 
exists a critical solidification velocity at which a bubble of 
given radius is engulfed. 161 Although these experiments have 
provided real-time observations of solidification dynamics, 
they do not necessarily reflect the behavior in a metallic sys- 
tem. Unlike salol, most metallic systems are nonfaceted. In 
addition, metallic systems have much higher interfacial energy 
compared to organic analogues. Further, the chemistry between 
a metal-gas system such as Al/H 2 is very different from an 
organic-organic system such as succinonitrile and acetone, 
where acetone is the source of pore formation. It must also 
be emphasized that experiments performed with organic ana- 
logues are typically carried out in thin slabs with character- 
istic thickness in the range of 10 to 300 pm. Consequently, 
the effect of wall friction on pore nucleation and distribution 
can be significant, resulting in controversial interpretations. 

The experimental objective of this study is to isolate the 
interaction between a gas pore and a planar SL interface in 
metallic systems. The understanding obtained from such an 
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investigation can be directly transferred to the distribution of 
interdendritic gas porosity and solutal segregation associated 
with it in castings solidification. To this goal, horizontal direc- 
tional solidification (HDS) experiments with pure A1 and Al- 
0.25 wt pet Au were conducted under direct and real-time x-ray 
transmission microscope (XTM) observation. This technique 
permits real-time tracking of the change in shape and size of 
porosity and the localized porosity-induced solutal segregation. 

Numerous models have been proposed over the years to 
model microporosity in castings. Piwonka and Fleming 171 
have elegantly modeled the interplay between gas and shrink- 
age porosity. They essentially showed that the effect of 
dissolved gas and shrinkage on pore formation is additive. 
Higher surface tension of the melt, on the other hand, can 
significantly increase the critical amount of dissolved gas 
required for gas porosity formation. Poirier et a/J 81 improved 
the model by calculating the local pressure in the mushy 
zone as a function of columnar dendritic microstructural fea- 
tures and the corresponding permeability and fluid flow in 
the interdendritic region. Lee and Hunt [4,9] used a continuum 
model to describe their experimental results in the Al-Cu 
alloy. They allowed pore growth by consumption of hydro- 
gen from a supersaturated control volume. At each time- 
step, the pore radius is compared to dendritic spacing to 
determine the shape of the pore. Voller and Sundarraj 1101 and 
later Rousset et al) n] modeled the effect of microporosity 
on solutal segregation by coupling an inverse segregation 
model with a microporosity model. They demonstrated that 
since inverse segregation is directly related to interdendritic 
fluid flow, the presence of microporosity decreases the extent 
of inverse segregation in the mushy zone. However, none 
of these models describes the exact physics that correlate 
solidification velocity with the SL interface shape as it inter- 
acts with a pore and the localized solutal segregation induced 
by the pore. In this article, we will address these issues from 
a modeling standpoint. 

The theoretical problem is to accurately track the sharp 
and arbitrarily shaped SL interfaces. Finding an adequate 
solution to this problem is an important and necessary step, 
because it can open new possibilities to study the microscale 
phenomena occurring during the solidification. Morpho- 
logical stability, coarsening processes, and SL interface inter- 
action with gaseous or non-metallic inclusions are only a 
few examples of phenomena that existing numerical solidi- 
fication models can describe with various degrees of accu- 
racy. Among the various techniques developed to date, the 
boundary-fitted coordinates phase field and 

level set 1 18,1 91 methods are the most notable. However, the 
boundary-fitted coordinates method experiences difficulties 
when dealing with highly distorted interfaces. It requires 
remeshing as the interface changes its shape. The versatility 
of the phase-field method lies in the fact that it can use a 
simple finite-difference technique, use a fixed grid, and does 
not require remeshing. However, despite tremendous efforts 
to improve the computational efficiency, 120,211 the method 
is still computationally intensive and produces only diffuse 
interfaces. The level set method was first introduced by 
Osher and Sethian 1181 and, since then, has been used to solve 
numerous problems 122-271 including solidification. The method 
is built on the idea that a moving interface separating two 
phases can be represented as the zero contour of a function 
a), which is the level set function. In general, the function 


c o is defined as the signed normal distance from the inter- 
face to each of the grid points on the computational lattice. 
It has its own equation of motion represented by 

da) 

— + F|Vw|=0 [1] 

dt 

where t is time and F is the normal velocity field, which is 
variable in space and time. A number of ways were pro- 
posed to construct the normal velocity field, F. A boundary 
integral approach was initially used by Sethian and Strain. 1191 
Then, faster algorithms based on partial differential equa- 
tions were developed. 1261 More recently, Kim et al [27] applied 
an even simpler method by calculating the velocity field 
from the difference in the normal derivatives of the tempera- 
ture field at the interface. In essence, the level set method is 
conceptually similar to the phase field method in that neither 
one explicitly tracks the interface. 

Udaykumar and co-workers [28, 29,301 have proposed an inter- 
face tracking technique based on a finite volume formulation. 
The moving SL interface is represented by a discrete set of 
markers that are positioned at the intersection of the interface 
with the grid lines. In this manner, after each iteration, only 
the interface must be reconstructed. However, the requirement 
of redefinition of the control- volumes near the interface con- 
siderably increases the difficulty of model implementation 
and makes it susceptible to numerical errors. 

In the present article, we propose a numerical technique 
based on the Udaykumar and Shyy model. It uses the finite- 
difference formulation, thus eliminating the need for control- 
volume redefinition at the interfaced 311 The proposed approach 
requires calculations of the normal velocity only at the SL 
interface. This is in contrast to the level set method where 
extension and reinitialization of the normal velocity off the 
interface is necessary. 

II. MODEL FORMULATION 

A. Interface Tracking 

Consider a SL interface, as shown in Figure 1. At the 
beginning of a time-step, the interface is described by a 
discrete series of markers positioned at the intersection of 
the SL interface with the grid lines. The open circles in Fig- 
ure 1 represent these markers. During the time- step, the SL 
interface advances according to the relative values of the 



j A new interfacial points D new interfacial points 
j- J}?i ojjjhe grid lines _ _ positioned on the gridjines_ 


Fig. 1 — Definition of the interfacial points. 
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temperature and solute gradients in the liquid and solid 
phases. These gradients and the corresponding local solidi- 
fication velocity are calculated at the marker points by fol- 
lowing a procedure that will be explained latter. At the end 
of this new time-step, the SL interface markers will find 
themselves in positions other than the intersection with the 
grid lines. These are represented by the open triangles in 
Figure 1. Consequently, a new set of markers, on the grid 
lines (the open squares), must be generated based on the 
positions of the old marker points. This task is done by a 
simple intersection procedure. This procedure requires the 
marker points to be numbered in increasing order, such 
that the liquid lies always to the left and the solid to the 
right as one advances along the interface in the direction of 
the arrow, as shown in Figure 2. Then, to calculate the inter- 
section of the SL interface with the grid lines passing between 
the points i-l and /, a variety of methods can be employed. 
Udaykumar and Shyy [30] have used segments of circular arcs 
to describe the SL interface between three consecutive 
markers. In a previous work by two of the authors of this 
paper, [32] second-degree polynomials were fitted through 
three consecutive markers, while in the present formulation, 
the cubic spline interpolation procedure [33] is used to find 
the position of the new marker points. 


B. Local Interfacial Gradients/Fluxes and Local 
Solidification Velocity 

The local solidification velocity, Le., the local rate at which 
the SL interface moves, is determined by the interfacial 
fluxes of the transport variables (heat and solute) calcu- 
lated at the marker points. Since these fluxes are normal to 
the SL interface, one needs to calculate first the normal direc- 
tion at each marker point. The direction of the unit normal 
vector, N, to a curve y = y(x\ is given by 


V 1 4- (dyldx) 2 



+ l*j 


[ 2 ] 


where the subscript i denotes that the calculation is perfor- 
med at the interfacial point j, and i and j are the unit vectors 
along the .r and y directions, respectively. Also, the function 
y(jc) can be constructed as a cubic polynomial fitted through 
four consecutive interfacial marker points. Because the inter- 
facial marker points are numbered according to the rule 
established in the previous section, the normal calculated 
through Eq. [2] will point always into the liquid phase. By 
simply changing the sign in Eq. [2], the normal will point 
into the solid phase. 

Consider now an interface segment containing the marker i. 
The gradient of a field variable (temperature, solute con- 


liquid 



centration) at point i is calculated based on the following 
equation (also Figure 3): 

^, + A W - <Ej R1 

dn An J 

where An is the length of a probe projected from the inter- 
facial point i along the direction of the normal vector N into 
the adjacent liquid and solid phases, and and 316 
the values of the field variables at the locations i and i 4- An , 
respectively. The problem when using Eq. [3] is that the 
location i 4 An usually does not fall onto a grid point, and 
therefore, the value must be obtained by using some 
kind of interpolation technique. The interpolation function 
for plane bilinear isoparametric elements, very popular in the 
finite-element method (for example. Reference 34), has been 
used in this work. This type of interpolation function is of 
the form <I>(x,v) = a 4- b*x 4 c-y 4 d-x-y, where the coeffi- 
cients a, b , c, and d can be obtained by choosing four 
grid points in the same phase around the location i 4 An 
(Figure 4). Then, the field variable at the location i 4 An 
characterized by the local coordinates (£ m 7j„) is given by 

4 

f 4] 

p= i 



Fig. 3 — Calculation of the interface gradient of a field variable. 



Fig. 4 — Development of the interpolation function. 
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where <& p is the value of the field variable at the node p 
( p = 1-4, Figure 4), and N p represents the shape functions 
defined as [32,34J 

' n , = i(l - £ n )(l - T ,„) 

Ni = j(l + f»)(l - T»«) 

< [5] 

V 3 = ^(1 - f«Xl + Vn) 

,N 4 = + f«Kl + Vn) 


A necessary requirement for interpolation is that the four 
grid points surrounding the location i 4- An must pertain to 
the same phase. In order to ensure this requirement, the 
length, An, of the projected normal is chosen as 

An = VA* 2 + Ay 2 [6] 


where Ax and Ay are the grid size along the x and y direc- 
tions, respectively. 

The interfacial temperature, T h at a location i is given by 
the relationship 

Ti =T M + m L > (Cdi + T • K t [7] 


where T M is the melting temperature of the pure base metal, 
m L is the liquidus slope, (C L ), is the solute concentration at 
the interface in the liquid phase, T is the Gibbs-Thomson 
coefficient, and K t is the SL interface curvature at the location /. 
The term K ( can be calculated as 


(d 2 y/dx 2 )j 
[1 + (dy/dxfi ] 3/2 


[8] 


where, again, y = y(x) is the local shape of the SL inter- 
face. The interfacial solute concentration in the solid phase, 
(C s )n is related to the liquid concentration, ( C L ) h through 
the relationship 


(C L ), = 


(Csh 

k 


[9] 


where k is the solute partition coefficient between the solid 
and the liquid phase. 

The local solidification velocity, (V n ) h normal to the inter- 
face at each marker point, /, can be calculated from heat and 
mass balance at the interface, as follows: 


{V n )iC p {dT) t = k L 



+ (V„)/A H f +k s 



[ 10 ] 


-d l • = (V n ), (1 - k) ■ (C L ), [11] 

where C p is the volumetric specific heat of the solidifying 
alloy; k L and k s are the thermal conductivities in the liquid 
and solid phase, respectively; A/// is the volumetric heat of 
fusion; D L is the solute diffusion coefficient in the liquid; 
and (d<bj/dri)i is the interface gradient of the field variable 
(temperature, T, and concentration, Q in the phase j (solid, S, 


and liquid, L). Note that the sign convention used for the 
gradient terms in Eq. [10] is defined by Eq. [3]. Also, in 
Eq. [10], the quantity {dT)i represents the change of the 
interface temperature at the location i in a time dt as a result 
of the interface curvature and concentration change during 
the time dt. The effect of solute diffusion in the solid phase 
has been ignored in Eq. [11]. 


C. Calculation of Field Variables 

The field variables can be calculated by solving the trans- 
port equations separately in the liquid and solid phases. Thus, 


dt 


H V 2 d> 


[ 12 ] 


where H is the generalized diffusivity, i.e., thermal diffu- 
sivity in the solid, a 5 , and in the liquid, a L , as well as mass 
diffusivity, D s and D L , in the solid and liquid, respectively. 
Equation [12] can be solved by using a finite difference 
method. 

When discretizing Eq. [12], a nonuniform grid size must 
be considered. Even if a uniform grid is used to discretize 
the computational domain, a nonuniform spacing appears in 
the proximity of the SL interface. This situation is illustrated 
in Figure 5. By using an implicit formulation, the finite 
differences expression of Eq. [12] can be written as 

= a • [b w < pf-u + M . + bt&l+tj 

+ - (b x + b y Wj [13] 

where the superscripts n and o indicate the value of the field 
variable at the end and at the beginning of the time-step, 
respectively, and A t is the size of the time-step. With the 



Fig. 5 — Nonuniform spacing in the vicinity of the SL interface. 
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notations used in Figure 5, the coefficients b w , b N , b& b s , 
b x , and b y in Eq. [13] are as follows: 


b w 

&E 

bs 

b N 

b x 


2Ax E IAx w 

Ax £ Ax^(1 + Ax E l Ax w ) 

2 

Ax £ Ax w (l + Ax £ /Ax w ) 
2 Av v / Ay$ 

A>VvA;ys(l + A\'aV A y^) 

2 

A^A>’ 5 (1 + Ay^/Ay*) 


Ax £ Axw 3 Ay N Ay s 


[14] 


To solve Eq. [13], boundary conditions on the SL inter- 
face and sample walls are required. Here, we must point 
out that Eq. [7] through Eq. [9] do not unequivocally define 
the boundary conditions on the interface, unless the solidi- 
fication of a pure substance is considered. The reason is 
that these equations do not account for the change of the 
interface concentration, ( C L ) h during the solidification 
process. Except for simple cases, such as directional solidi- 
fication with a planar SL interface. [35] there are no ana- 
lytical solutions to describe this change. However, to 
account for the evolution of ( C L ) h the following method- 
ology has been used in the present work. It relies on the 
assumption that the local quantities O l+An (temperature, 
concentration) as well as the local curvature, K n do not 
change significantly as a segment of the SL interface 
moves over the distance (V„)rAt. Therefore, the local 
change of the SL interface temperature, A T h during the 
time At is given by 


A T, = m L [(C L )i - (cy?j [15] 

where (C £ )° is the concentration at the beginning of the 
time-step. In addition, the finite difference expressions 
for the balance equations (Eq. [10] and Eq. [11]) can be 
written as 


(V n )i C P AT, 


h[r nL 


(T? + fATd\ 
An 


(V„), A H f 


ksttTf + fATd - T aS ] 
An 


[16] 


(Vn)i = D l 


( Cdi ~ c nL 

An(l - k)(C L )i 


[17] 


where T ,° is the local interface temperature at the beginning 
of the time-step; T nL and C nI are the temperature and con- 
centration, respectively, at the location i + An into the liq- 
uid phase; T nS is the temperature similarly defined in the 
solid phase; and / is a factor that takes the value of 0.5 to 
account for average temperature gradients at the interface 
during the time-step. Further, by eliminating (CJ, from 
Eq. [17] and replacing it in Eq. [16] by making use of 
Eq. [15], the following second degree equation can be obtain 
for the local growth velocity (V,,),-: 


MVJ - A 2 (V n ), + A, = 0 [18] 


with 


r 

A\ 

A-2 

< a 3 

Q 


= l A H f + Cpm L (C L )° J(1 - k)An 
= [A H f + CpmMCtX - C nL )]D L - Q( 1 - k)An 
fm L {k s + k L )C nL 

An " Q \° l 

= f>n, (k s + k L )(C L )° T nL - T[ _ T° - T nS 
An L An 5 An 


[19] 


The only valid solution of Eq. [18] is 


(V«)i 


Aj A 2 AA\A 3 

2A X 


[ 20 ] 


Then, (C L ), and T t can be calculated from Eqs. [17] and [7], 
respectively. The boundary conditions on the SL interface 
and the solidification velocity are thus calculated. It must 
be added here that although the probability is very small, 
in some cases, it is still possible for an interface marker 
point to fall on a grid point, thus making at least one coeffi- 
cient in Eq. [14] take an infinite value. For such special 
situations, the particular grid point located on the SL interface 
is treated as a boundary point. While the temperature of this 
boundary point is the same on either side, the concentration 
is different on the side facing the liquid phase from that on 
the side facing the solid phase (Eq. [9]). 


D. Nondimensional Equations and Time-Step 

In order to minimize the round-off errors, a nondimensional 
form of the constitutive equations presented in this section 
must be used. Because for all the case studies considered in 
this work an initial solidification velocity, V, was imposed 
through the external boundary conditions, this velocity was 
taken as a reference velocity, V ref . The reference time scale, 
tref, then follows from the relationship f ref = where L^ f 

is the characteristic lengthscale of a particular case. The other 
reference quantities used in this work are the melting/solidus 
temperature of the solidifying material, liquid thermal and solu- 
tal diffusivities, liquid thermal conductivity, and the initial 
concentration of the melt. 

Because the SL interface is explicitly tracked in this 
model, restrictions must be imposed on the size of the time- 
step in order to ensure stability. A variable time-step of size 
given by the following empirical relationship has been used 
in this work: 


At = 


2 ‘ Vn, 


Ay 3 • V 


[ 21 ] 


where F ITiax is the instantaneous maximum solidification veloc- 
ity in the system and V is the imposed solidification velocity. 


III. MODEL INTERROGATION 

In Sections A and B, a number of case studies are considered 
in order to test the capabilities and the accuracy of the proposed 
numerical model. Whenever possible, analytical solutions will 
also be presented for comparison. The thermophysical proper- 
ties used in the calculations are presented in Table I. 
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Table I. Thermophysical Data Used in Calculations 


Property 


A1 Al-2 Wt Pet Cu 

Al-0.25 Wt Pet Au 

Sn-0.5 At. Pet Bi 

i Zr0 2 

H 2 

Heat of fusion 141,421 (J/m 3 )*10 8 


9.5 

9.5* 

9.5* 

4.179* 

— 

— 

Gibbs-Thomson coefficient 1411 








(K-m)-lCT 7 


0.9 

0.9* 

0.9* 

— 

— 

— 

Liquidus slope [4 1,43,441 (K/wt%) 


— 

-2.6 

-1.43 

-1.95** 

— 

— 

Partition coefficient 141,431 


— 

0.14 

0.048 

0.28 

— 

— 

Solute diffusivity in the 








liquid 141,451 (m 2 /s)T0 -9 


— 

3 

8 

1.1 

— 

— 

Thermal conductivity 141,42,46,471 

liquid 

95 

95* 

95* 

30 

1.97 

0.43 

(W/m-K) 

solid 

210 

210* 

210* 

60 



Specific heat 141,42,47,481 

liquid 

2.58 

2.58* 

2.58* 

1.81* 

3.63 HO 6 

387.92 

(J/m 3 *K)T0 6 

solid 

3 

3* 

3* 

1.96* 



♦Property assumed to be the same as that of the base metal. 
**K/at. pet. 


A. Solidification of Pure Aluminum in the Vicinity 
of a Particle 

Directional solidification of pure aluminum in the vicinity 
of a spherical Zr0 2 particle of radius R p = 22.5 ^m was 
chosen as a first case study. The imposed solidification 
velocity was V = 2 /zm/s under a thermal gradient in the 
liquid G l = 1 K/mm. The calculations were performed on 
a sample of length L — 2 mm and a width W, which was 
varied, from case to case, from 1 to 0.4 mm in order to test 
its influence on the evolution of the SL interface shape. 
However, no noticeable effect was observed. At the begin- 
ning of the calculations (Le, y t = 0), it was considered that 
at one end of the sample a slice of Az = 10 fim thickness 
was in solid state and that the SL interface was perfectly 
planar. The center of the particle was placed at x = Wf 2 and 
z = Az + h cy where h c is the distance between the particle 
center and the planar SL interface (Figure 6). Insulated 
boundary conditions were used on the x direction, while at 
the two ends of the sample, constant cooling rates were con- 
sidered (also Figure 6). The cooling conditions were set such 
as to maintain the imposed solidification velocity, V, of the 
originally planar interface for the entire duration of the cal- 
culations. They were dTIdt — ~V G L at the liquid end and 
dTIdt = — V G S at the solid end, where G s is the tempera- 
ture gradient in the solid phase. 

The available analytical solution for a problem of this 
type [36] is presented in Appendix I together with some scaling 
considerations. This solution predicts that when the thermal 
conductivity of the particle, k Py is smaller than that of the 
liquid, k Ly the SL interface will become convex toward the 
particle, forming a bump as it approaches the particle. Indeed, 
our numerical calculations show this behavior, as can clearly 
be observed in Figures 7 and 8. 

Figure 7 shows a comparison between the numerical and 
the analytical solutions when neglecting the Gibbs-Thomson 
effect (i.e., T = 0). The agreement is good, thus proving the 
capability of the proposed interface tracking model. In all 
our numerical simulations, the interface perturbation is some- 
what sharper than that predicted by the analytical solution. 
This is mainly because of the insulated boundary conditions 
used at a finite distance from the particle (Figure 6) as com- 
pared to the analytical solution, which implies that dT/dx = 0 
occurs at x — > oo. 
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Fig. 6 — The geometry and boundary conditions for the sample used in 
computations. 

Figure 8 shows the shape of the SL interface when account- 
ing for the Gibbs-Thomson effect. It can be observed that, as 
expected, the amplitude of the interface perturbation becomes 
considerably smaller as compared to the case when T = 0. 
The influence of the temperature gradient, G L , is also shown 
in this figure. In qualitative agreement with Eq. [Ic] (Appen- 
dix I), the numerical calculations show that, for the same 
distance h c , as G L increases, the perturbation becomes sharper 
and its amplitude also increases. This implies that the pertur- 
bation develops faster when increasing G L . 

It also has to be mentioned here that when accounting for 
the Gibbs-Thomson effect, the outcome of the calculations 
can be influenced to some extent by the choice of the grid 
size on which the calculated curvature is strongly depen- 
dent. A procedure to estimate the level of errors introduced 
by a certain grid size is presented in Appendix II. Although 
the developed model allows the use of a nonuniform grid, 
a uniform grid of size R P l 6 on both x and y directions was 
used to produce the preceding results. For all the other prob- 
lems that will be discussed in this article. This paper unless 
specified otherwise, the same grid size has been employed. 

B. Solidification of Alloys; the Influence of Solute 

The first test problem that we consider in this case is the 
solute redistribution during the initial transient of plane front 
directional solidification of a binary alloy. An analytical 
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Fig. 7 — SL interface shape: comparison of analytical and numerical solutions (pure Al, Zr0 2 particle, R p = 22.5 /xm, V ~ 2 /xm/s, with no Gibbs-Thomson effect). 
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solution to this problem 1351 shows that the concentration pro- 
file, C 5 (z), in the solidified sample can be expressed as 


Cs(z) 
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[ 22 ] 


where C 0 is the initial concentration in the liquid sample, z 
is the distance measured from the beginning of the sample 
along the solidification direction, and erf and erfc are the 
error and complementary error functions, respectively. 

The results shown in Figure 9 are for Sn-0.5 at. pet Bi 
alloy at V = 4.9 /xm/s and a thermal gradient in the liquid 
G L = 1 7.3 K/mm. The numerical calculations were performed 
for a sample of length L = 10 mm and width W = 20 /xm 
using a grid size of 4 /xm. In order to maintain the imposed 
solidification velocity and temperature gradients for the entire 
duration of the calculations, the following boundary condi- 
tions for the temperature at the two ends of the sample have 
been used: 

At the liquid end. 



Solidified distance, m 

Fig. 9 — Numerical predictions vs analytical solution 1351 for Bi segregation 
during the initial transient plane front solidification (V = 4.9 jxm/s) of 
Sn-0.5 at. pet Bi. 
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where dC[/dt represents the time change of the liquid con- 
centration at the SL interface. Insulated boundary conditions 
for both temperature and concentration were used on the x 
direction (Figure 6). 

A good agreement between the analytical solution and the 
numerical results can be observed in Figure 9. In general, 
the numerical solution gives slightly higher values for the 
concentration as compared to the analytical solution. This 
is because of the linear approximation used to calculate the 
concentration gradient at the SL interface. 

With our numerical model thus validated against analy- 
tical solutions and experimental results for the solidification 
of binary alloys, we can proceed to study the evolution of 
the SL interface shape when a foreign particle is added to 
the melt in front of the interface. We chose Al-2 wt pet Cu 
for the alloy and Zr0 2 for the particle. A steady-state distri- 
bution for the solute and linear distribution for the tempera- 
ture were considered as initial conditions. The boundary 
conditions for temperature are presented in Figure 6. Insulated 
boundary conditions on the x direction as well as on the 
surface of the particle were added for the solute. The results 
are shown in Figure 10 for three values of the distance, h c , 
between the center of the particle and the flat SL interface. 
It can be observed that in the initial stage of solidification, 
when the interface is relatively far from the particle, the inter- 
face forms a small bump because the thermal conductivity 


of the particle is smaller than that of the melt. At this stage, 
the solute diffusion is only slightly obstructed by the presence 
of the particle, and therefore, the outcome is mainly deter- 
mined by the thermal field. As the solidification proceeds, 
however, the solute diffusion becomes more and more affected 
by the presence of the particle. The localized solute accu- 
mulation decreases the interface temperature and the inter- 
face begins to turn concave beneath the particle. The depth 
of this concavity increases as the particle/interface distance 
decreases. Also note that for a very small advance of the 
flat SL interface, the amplitude of the two humps that form 
on both sides of the vertical centerline increases dramatically. 
This is because the lateral diffusion gradually intensifies 
in these regions, and therefore, the solute depletion occur- 
ring around the tip of the two humps allows them to grow 
ever faster. 

While for the solidification of pure substances our numer- 
ical calculations as well as the analytical solution 1361 show 
that the solidification velocity has a negligible effect on the 
evolution of the SL interface near a foreign particle, this 
conclusion cannot be extended to the solidification of alloys. 
The influence of the solidification velocity is shown in 
Figure 11. The solute buildup in the particle/interface gap 
becomes more significant with the increase of the solidifica- 
tion velocity, and therefore, the depth of the concavity formed 
on the SL interface beneath the particle will be higher at 



Fig. 10 — Time evolution of an alloy SL interface approaching a foreign particle (Al-2 wt pet Cu, Zr0 2 particle, R P = 22.5 jum, V = 0.2 fx m/s, G L = 7 K/mm). 
For a better view of the interface deformation the particle shape appears noncircular because of nonequal scales used on the vertical and horizontal directions. 
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Fig. 1 1 — Influence of the solidification velocity on the shape of the SL interface when approaching a foreign particle (Al-2 wt pet Cu, Zr0 2 particle, R P = 
22.5 ixm, G L = 7 K/mm). 
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higher solidification velocity. On the other hand, when meas- 
ured with respect to the flat region of the interface, the ampli- 
tude of the two humps that form on both sides of the vertical 
centerline is smaller at higher velocities. Thus, the solute 
effect becomes predominant over the thermal effect (quan- 
tified through the ratio k P lk L ) when increasing the solidifica- 
tion velocity. Note that for the velocity V = 2 jiim/s. used 
to generate Figure 11, the stability theory 1 ' 71 predicts that 
the SL interface is inherently unstable and perturbations 
other than those induced by the presence of the particle 
should also appear on the interface. This is true and we have 
demonstrated 1381 that the present model has the capability 
to describe the development of such perturbations. For the 
case presented in Figure 11. however, where the sample 
width. W. was considered as only 0.4 mm. the perturbations 
induced by the presence of the particle override the pertur- 
bations that would have been induced by the constitutional 
undercooling effect alone. 

So far, our numerical model was validated against simple 
problems for w hich analytical solutions are available. It was 
also show n that for more complicated problems, the model 
provides meaningful results. However, a complete validation 
can be performed only against accurate experimental data. 
In Section IV, we will present such an experiment and 
compare the numerical predictions with the experimental 
measurements. 


IV. EXPERIMENTAL APPROACH 

The principle of the experimental setup is presented in 
Figure 12. Projection radiography using a microfocus X-ray 
source offers magnification, adequate resolution, and. with 
suitable detector technology, adequate contrast. A state-of- 
the-art submicron source w ith acceleration voltage of 10 to 
100 kV allows imaging of solidification of metal alloys. 
Real-time viewing is possible by employing an X-ray image 



Fig. 12 — A schematic illustration of the XTM apparatus and principle. 140 


intensifier and a cooled CCD camera. Resolution, which is 
limited bv the X-ray spot size, can approach micrometer 
values. Radiography by projection permits placing the speci- 
men in a furnace between the X-ray camera and the source. 

The XTM furnace employed in this research is a modified 
horizontal Bridgman-Stockbarger design operated in air or 
inert gas. Samples with typical dimensions of 19 X 4 X 
1 mm are placed in the cavity of an X-ray transparent boron 
nitride crucible. The sample is then covered with a boron 
nitride lid and placed in the furnace. A stepper motor driven 
screw translator translates the specimen in its crucible at 
rates of 0.1 to 100 gm/s. During solidification, using the 
XTM. it is possible to track the interface position as a func- 
tion of time and thereby obtain true interface velocity, as 
opposed to furnace velocity. 

For this investigation, the selected sample materials were 
pure A1 and Al-0.25 wt pet Au. The Al-Au alloy was selected 
because the significant density difference between the solute 
and the solvent provides good X-ray contrast. The samples 
were melted in an induction furnace and cast into cylindrical 
rods. LECO* combustion measurements of the as-cast samples 

*LECO is a trademark of LECO Corporation. St. Joseph. MI. 

showed a H : content of 0.29 ppm in the pure A1 sample and 
0.3! ppm in the Al-Au sample. The cylindrical castings were 
then cold rolled into 1-mm-thick strips and trimmed to fit 
into the cavity of the boron nitride crucible. The HDS exper- 
iments were performed to investigate the dynamics of the 
interaction between a pore and the SL interface. The 
processed samples were sectioned and the solute segregation 
pattern in the vicinity of a pore w as quantified using micro- 
probe analysis. 

A. Pure Aluminum 

The in-situ radiographic observation of pore growth in 
liquid A1 is shown in Figure 13. The translation velocity 
w as set at 4 jiim/s and the thermal gradient ahead of the inter- 
face was 4.7 K/mm. Prior to its interaction w ith the SL inter- 
face. the pore was observed to grow w ith a circular cross 
section (Ligure 13(a)). When the SL interface makes contact 
with the pore, the pore grows as an ellipsoid of revolution 
with its major axis aligned with the solidification direction 



(a) (b) (c) 


Fig. 13 — Real-time radiographs of a pore and of the SL interface moving 
leftward: (a) prior to interaction with the SL interface: and (h) and (c) sub- 
sequent ellipsoidal evolution during interaction with the SL interface in 
pure Al. The white arrow points to the position of the SL interface. 
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(Figures 13(b) and (c)). Figure 14 shows the measured time 
change of the volume and aspect ratio of a pore. The dashed 
vertical line indicates the time at which the interaction 
between the pore and the SL interface is initiated. It is evi- 
dent from Figure 14 that the growth rate of the pore after it 
starts interacting with the interface is significantly higher 
compared to the growth rate prior to this interaction. Similar 
growth behavior of gas pores has also been observed in trans- 
parent organic alloys. 151 The sudden increase of the growth 
rate of the pore implies that more hydrogen is available for 
its growth than that available through bulk diffusion in the 
melt. The pore acts as a sink for the H 2 rejected at the solidifi- 
cation front. 

B. Al-0.25 wt pet Au Alloy 

Figure 1 5 shows a postsolidification radiograph of a pore 
after engulfment by a planar SL interface. The pore has an 
ellipsoidal shape with a size of about 800 gm along the 
major axis and 660 pm along the minor axis. The growth 
velocity was 5 pm/s and the temperature gradient 4.7 K/mm. 
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Fig. 14 — Experimentally measured change in pore volume and aspect ratio 
as a function of time in pure Al. 


The purpose of this experiment was to investigate the effect 
of the pore on the solutal segregation pattern. The radiograph 
shows that there is a strong segregation of Au around the 
pore just prior to engulfment. The scanning electron micro- 
graph (Figure 15(b)) shows that the microstructure in this 
area consists of eutectic colonies. The eutectic structure 
inside one such colony is shown in Figure 15(c). The thick- 
ness of this eutectic “band” varies between 70 and 80 pm. 
The centerline microprobe scan indicates that the Au com- 
position in this region is approximately 8 wt pet. This is in 
agreement with the equilibrium phase diagram for the Al-Au 
system, which shows a eutectic transformation at 7.46 wt pet 
Au. Subsequent to engulfment of the pore, a comet tail- 
shaped segregation pattern is clearly evident on the left of 
the pore (Figure 15(a)). The mechanism of formation of this 
strongly segregated region will be explained in Section V. 


V. NUMERICAL CALCULATIONS 
VS EXPERIMENTAL MEASUREMENTS 

At the outset, it should be mentioned that at the present 
stage, the model does not have the capability to account for 
changes of pore size and shape during the interaction with 
the SL interface. Therefore, the pores have been assumed 
to be spherical and of constant radius. 

Figure 16 shows the calculated shape of the SL interface 
in various stages of a H 2 pore engulfment by the solid phase. 
The occurrence of a trough on the SL interface after engulf- 
ment is observed. This is typical, as suggested by Eq. [Id], 
for the situations when the thermal conductivity of the par- 
ticle (H 2 pore in this particular case) is smaller than that 
of the metallic matrix. As can be seen in Figure 16, once 
the pore is completely engulfed, the two segments of the 
interface are again joined together. The model handles 
this special situation as follows: on either side of the pore, 
the coordinates of the intersection of the SL interface and 
the pore are tracked during the engulfment of the pore. In 



Fig. 15 — (a) through (c) Postsolidification radiograph showing segregation of Au around the porosity. Superimposed is a centerline microprobe scan. 
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Distance along the sample width, m 

Fig. 16 — Numerically calculated shape of the SL interface in various stages 
of engulfing a H : pore of fixed size (pure Al. = 22.5 ^m. V 7 = 4 fi m/s, 
G l = 4.7 K7mm). 
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Distance along the solidification direction, m 

Fig. 17 — Numerically calculated Au distribution profile along the sam- 
ple centerline showing solute trapping in the SL interface/pore gap 
(Al-0.25 wt pet Au. V 7 = 5 /jliu/s. G L = 4.7 K/mm. R^ Te = 330 /xm). 

this model, the adopted criterion for the complete engulf- 
ment is the crossover of these two points. When the 
crossover occurs, a new interface marker is generated at 
the point of intersection. 

For the solidification of the Al-0.25 wt pet Au alloy, our 
calculations also predict that at a distance from the pore 
similar to that observed experimentally, the Au concentra- 
tion raises to the eutectic value, as shown in Figure 17. At 
the time of this publication, a model for eutectic solidifica- 
tion was not yet implemented in the numerical program, and 
therefore, the computation was stopped when the Au con- 
centration in the liquid pore/SL interface gap reached the 
eutectic value. This precluded a quantitative prediction of 
the second region of strong Au segregation developed in 
front of the pore after its complete engulfment. Neverthe- 
less, based on the interface shape of pure Ai (Figure 16). it 
is reasonable to assume that prior to the complete engulf- 
ment. a dramatic increase of solute concentration occurs in 
front of the pore. This is because the solute diffusion is now 
directed toward the centerline of the sample, as shown in 
Figure 18(a). Therefore, the solidification in this region is 
delayed. Consequently, the shape of the SL interface at the 
moment of complete engulfment will look as show n in Fig- 



Fig. 18 — The mechanism of the comet tail-shaped segregation region in 
front of the H : pore: (a) solute diffusion toward the sample centerline before 
complete engulfment: and ( b ) solute trapping in front of the pore at the 
time of complete engulfment. 


ure 18(b). thus explaining the comet tail-shaped region of 
high Au concentration, which appears in Figure 15. 

VI. SUMMARY 

Real-time experimental observations of the interaction 
between a FL pore and a planar SL interface as well as a 
numerical mathematical model that explains these observa- 
tions have been presented in this article. 

The real-time measurements of pore growth in pure Al have 
revealed the existence of two growth regimes, as follows. 

1 . When the pore is relatively far from the SL interface, the 
mechanism of pore growth is the hydrogen diffusion 
through the liquid phase. The growth rate increases as 
the SL interface approaches the pore. At this stage, the 
pore has a circular cross section. 

2. A sudden increase of the pore growth rate occurs when 
the solutal (i.e.. hydrogen) field ahead of the SL inter- 
face begins interacting with the pore. The pore begins 
to assume an ellipsoidal shape as the SL interface grad- 
ually engulfs it. 

For the solidification of an Al-0.25 wt pet Au alloy, both 
the X-ray imaging and the microprobe analysis have demon- 
strated the development of two distinct regions of high solute 
concentration around the pore. The two regions have dif- 
ferent shapes because the mechanism of their formation is 
also different. Chronologically, the first region develops in 
front of the pore as a result of solute built up in the pore/SL 
interface gap. It follows the contour of the pore and. there- 
fore, develops a moonlike shape. The second region devel- 
ops behind the pore, where, initially, the solidification is 
delayed mainly because of the thermal conductivity effect 
(i.e.. the thermal conductivity of the pore is smaller than that 
of the solid phase growing around it). The solute diffusing 
toward this region further accentuates the solidification delay. 
Thus, the SL interface is forced to detach from the pore 
w hile it continues to advance at a certain angle tow ard the 
major pore axis where bridging eventually occurs. The liq- 
uid pool left behind the location where the SL interface 
bridging takes place is rich in solute and has a conical shape. 
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Consequently, porosities in castings act not only as material 
discontinuities that negatively affect the mechanical proper- 
ties, but also contribute to the development of regions of 
strong solute segregation, which further accentuates the detri- 
mental effect on the casting properties. 

A numerical mathematical model aiming to explain the 
experimental observations has been developed. It has the 
capability to track sharp and arbitrarily shaped SL inter- 
faces during solidification of pure substances and binary 
alloys. Any local interface temperature change because of 
Gibbs-Thomson effect or solute redistribution between the 
solid and the liquid phase can be accommodated. For a num- 
ber of relatively simple problems, the numerical predictions 
proved to be in good agreement with the available analytical 
solutions. 

With regard to the interaction between the SL interface 
and a nonmetallic inclusion (gas pore or ceramic particle), 
the model can accurately describe the change of the interface 
shape during this interaction. It was shown that the local SL 
interface dynamics in the vicinity of a nonmetallic inclusion 
are influenced not only by the thermal conductivity ratio, 
kplk L , but also by the imposed temperature gradient and solidi- 
fication velocity. The solute segregation pattern developing 
around such an inclusion can also be described by the model. 

The numerical calculations were performed on a LINUX 
platform with AMD Athlon 2000 MP processor, 1 .67 GHz. 
For a typical grid size of Ay = 3.5 yum, V = 2 ^m/s, and a L = 
37- 1CT 6 nr/s, the average computation time, required to simu- 
late 1 seconds of real time, is approximately (1.8-/V- 10“ 3 + 6) 
seconds, where N is the total number of the grid points. 
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APPENDIX I 


Analytical solution for the shape of the solid/liquid 
interface in the vicinity of a particle 

Consider a liquid phase of thermal conductivity k L in which 
the thermal field is characterized by constant temperature 
gradient, G L . If a foreign particle of radius R P and thermal con- 
ductivity k P is added to this liquid, then, for a two-dimensional 
case, the thermal field around this particle can be described 
by the relationship 1361 


T(x, z) = T, 


M 


7 _ k . R 2 K-z 

4 k a 2 , n 

x + ( h c — 


zf 


[la] 


where T is temperature, jc and z are the cartesian coordinates, 
T m is the melting point, and k a is given by 


k 


a 


+ k P 


[lb] 


Also, in Eq. [Ia], it was considered that the center of the 
particle is located at the position x = 0 and z = h c (Eq. [la]). 
Therefore, in the absence of the particle (or if k P = k L \ the 
temperature in the plane situated at z = 0 is equal to Tm, 
that is, the SL interface coincides with this plane. 


For the solidification of a pure substance, the local temper- 
ature, Tj(x), of the SL interface is given by the relationship 

Tj(x) = T M + T - K(x) [Ic] 


where T is the Gibbs-Thomson coefficient and K(x) is the 
local curvature of the interface. From Eqs. [Ia] and [Ic], it 
follows that the shape, Z 7 (x), of the SL interface can be cal- 
culated with the equation 


Zi(x) - 


ka'Rp 


K - Zj(x) 

x 2 + [h c - Z/U)] 2 


r - K(x) 
Gl 


[Id] 


Equation [Id] predicts that if k P 9^ k L , the SL interface will 
be planar only at x — > while in the vicinity of the particle, 

it will be concave or convex, depending on the ratio k P /k L . 
If k P /k L < 1 (i.e., k a > 0), the interface will be convex toward 
the particle, as shown in Eq. [Ia], and concave if k P /k L > 1 
(/.£., k a < 0). It turns out that, because of the presence of 
the curvature term, K(x ), which is given by Eq. [8], the 
general solution of Eq. [Id] is difficult to find. However, for 
a sufficiently high temperature gradient, G L , the right-hand 
side of Eq. [Id] becomes negligible and an approximate 
solution can easily be obtained. Obviously, this solution is 
exact only for the imaginary case when T = 0 (i.e., no 
Gibbs-Thomson effect). 

Inspection of Eq. [Id] suggests that by choosing the refer- 
ence length scale 

Pcf = Vfcj ' R P [Ie] 

where \k a \ is the absolute value of the parameter k a , Eq. [Id] 
can be written in nondimensional form as 


m = 


G L \k a \Rp 


Z, (x) - sgn (kj ■ 


K - z,(jf) 

(x ) 2 + (h c -z,(x)) 2 _ 

[If] 


The sign “ A ” in Eq. [If] indicates nondimensional quantities, 
while sgn (k a ) means the algebraic sign of k a . 

Equation [If] shows that the influence of the material and 
process parameters (i.e., G L , k a , T, and R P ) on the local SL 
interface curvature can be completely separated from the local 
interface coordinates when choosing the reference length- 
scale, as given by Eq. [Ie]. Thus, it seems that the natural 
lengthscale for this type of problems is not particle radius, 
R P , but V&p. This reference lengthscale has been adopted 
in this work. 



Fig. Ia — Nonplanar SL interface in the presence of a foreign particle 
( k p < k L ). 
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APPENDIX II 

The effect of the grid size on the calculated SL 
interface curvature 

One of the drawbacks of the interface tracking method 
presented in this article is that the accuracy to which the local 
SL interface curvature is calculated depends on the choice 
of the grid size. An illustration of the expected level of errors 
introduced by a certain grid size is presented in this Appendix. 

In Appendix I, it was shown that at a very high-temperature 
gradient, G L , or in the hypothetical case of Y = 0, an exact 
analytical solution can be obtained for the SL interface by 
means of Eq. [Id]. The corresponding exact curvature can 
also be calculated at any point from the same equation. It 
can be shown that the minimum interface curvature radius, 
Rj min, occurs when the SL interface first touches the particle 
at jc = 0 (Eq. [la]) and is given by 


R 


Amin 


K ~ 1 

K 


Rr 

2 


[Ha] 


For the particular case of an H 2 pore in a liquid A1 matrix 
(k a = 0.991) and R P = 22.5 /xm, Eq. [Ila] gives a curvature 
radius Rj = 0.102 ^m. This is a very disadvantageous case, 
because in order to minimize the errors introduced by the 
numerical calculation of the curvature, a grid size smaller 
than 0.1 gm would be necessary. In reality, however, the 
Gibbs-Thomson effect is always present and the result is that 
Rp min is larger than that predicted by Eq. [Ila]. A recently 
developed procedure 1391 shows that for this particular case, 
a more realistic value is R I mi n — 8 fim, which closely 
corresponds to the hypothetical case of T = 0, k P = 
24.92 W/m*K, and keeping unchanged all the other para- 
meters. Thus, even for this very disadvantageous case (i.e., 
k P k L ), the minimum interface curvature radius is no 

less than Rpl 3. 

Further, in order to test the influence of the gird size on 
the curvature calculated by means of the cubic spline method, 
the following procedure has been used in this work. 

1. At each location, x = (X, 4- X /+] )/2, where X, is the posi- 
tion of the /th grid line; the interface position, Z /(JC) , is 
calculated from Eq. [Id], for T = 0 and the newly calcul- 
ated k P . Note that the grid size on the x direction is given 
by the difference (X /+1 — X f ), and the interface marker 
points thus calculated do not coincide with the grid lines. 

2. Use the cubic spline interpolation procedure to position 
the interface marker points on the grid lines and calculate 
the curvature, K h at each such point. 

3. Calculated the curvature error. Err,, as 


Err, = 100 • 



mb] 


where Kf is the exact curvature at the /th marker point, 
calculated from Eq. [Id]. 

An analysis of the errors thus calculated indicated that the 
highest deviations occured near the inflexion points of the 
interface and at the point of minimum curvature radius (/.<?., 
at x = 0 in Eq. [Ia]). The deviations at the inflexion points 
can be as high as 2.5 times larger than those occurring at 
the point of minimum curvature radius. However, near the 



Grid size, m 


Fig. Ha — Calculated curvature error at the point of minimum interface curva- 
ture radius (R,, = 22.5 junt, k L = 95 W/m-K, k p = 24.92 W/m-K, T = 0 Km). 

inflexion points, the curvature assumes relatively low val- 
ues, and therefore, the numerical errors introduced at these 
points do not significantly influence the interface motion. 
Consequently, the attention should be focused at the point 
of minimum curvature radius. A plot of the curvature error 
occurring at this point vs grid size is presented in Figure [Ha]. 
Thus, for a grid size smaller than Rp/6, i.e., smaller than half 
of the minimum curvature radius. Fig. [Ila] indicates that 
the errors will be smaller than 8.5 pet. 


LIST OF SYMBOLS 


A u A 2 , A 3 , Q quantities defined in Eqs. [18] through [20] 
to calculate the normal growth velocity 
b E , b N , b s , coefficients in the discrete form of the 

b x , by, b w transport equation of field variables 


c 

concentration 

C L , C$ 

liquid and solid concentrations, respectively 

Cp 

volumetric specific heat 

Du D s 

solute diffusion coefficient in the liquid and 
solid phases, respectively 

<D 

generalized field variable 

/ 

coefficient that takes the value of 0.5 

Gu G s 

temperature gradient in the liquid and solid, 
respectively 

AH, 

latent heat of fusion 

h c 

distance between the planar interface and the 
center of the particle 

k 

partition coefficient 

k/ . k$, kp 

thermal conductivities of the liquid, solid, 
and particle, respectively 

a u «s 

thermal diffusivities of the liquid and solid, 
respectively 

K 

solid/liquid interface curvature 

r 

Gibbs-Thomson coefficient 

N u N 2 , 

shape functions used for the interpolation of 

n 3 ,n 4 

the field variables 

N 

unit normal vector 

An 

length of the probe projected along the normal 
direction 

m L 

liquidus slope 


METALLURGICAL AND MATERIALS TRANSACTIONS A 


VOLUME 35A, MAY 2004—1537 




Rp particle radius 

T temperature 

Tm melting temperature 

T u T s temperature in the liquid and solid phases, respectively 
t, At time, time-step size 

V imposed growth velocity through external boundary 

conditions 

V n normal growth velocity 

x, z Cartesian coordinates 

S generalized diffusivity 

Suhscrivts/Superscripts 
L pertaining to the liquid phase 

S pertaining to the solid phase 

i at the solid/liquid interface 

o at the beginning of the time-step 

nL, nS on the normal direction in the liquid and solid 
phases, respectively 
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